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ON THE CONVERGENCE OF SECOND ORDER SPECTRA AND 
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Abstract. Let A be a self-adjoint operator acting on a Hilbert space. The notion 
of second order spectrum of A relative to a given finite-dimensional subspacc C lias 
been studied recently in connection with the phenomenon of spectral pollution in 
the Galcrkin method. We establish in this paper a general framework allowing us 
to determine how the second order spectrum encodes precise information about the 
multiplicity of the isolated eigenvalues of A. Our theoretical findings are supported 
by various numerical experiments on the computation of inclusions for eigenvalues 
of benchmark differential operators via finite element bases. 
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1. Introduction 



Let ^ be a self-adjoint operator acting on an infinite dimensional Hilbert space "H and 
let A be an isolated eigenvalue of A. For I C R let 

lx{A) = I dE^ 

where E^ is the spectral measure associated to A. The numerical estimation of A whenever 

(1) Trl(_^,A)(^) = Trl(;,,^)(A) = oo, 

constitutes a serious challenge in computational spectral theory. Indeed, it is well es- 
tablished that classical approaches, such as the Galerkin method, suffer from variational 
collapse under no further restrictions on the approximating space and therefore might lead 
to spectral pollution P [4l [5l [T3l [ TT l l25l [23l [^ . 

The notion of second order relative spectrum, originated from [T3] (see Definition IT] 
below), and has recently allowed the formulation of a general pollution- free strategy for 
eigenvalue computation. This was proposed in [261 121) and subsequently examined in 
[5] [3 111! I28| . Numerical implementations of the general principle very much preserve the 
spirit of the Galerkin method and have presently been tested on applications from Stokes 
systems [21], solid state physics [10], magnetohydrodynamics [25; and relativistic quantum 
mechanics [5]. 

In this paper we examine further the potential role of this pollution-free technique for 
robust computation of spectral inclusions. Our goal is two-folded. On the one hand, 
we establish various abstract properties of limit sets of second order relative spectra. 
On the other hand, we report on the outcomes of various numerical experiments. Both 
our theoretical and practical findings indicate that second order spectra provide reliable 
information about the multiplicity of any isolated eigenvalue of A. 

Section [2] is devoted to reformulating some of the concepts from [261 1211 [G] [7] [11] 
allowing a more general setting. In this framework, we consider the natural notion of 
algebraic and geometric multiplicity of second order spectral points (Definition [2| and 
establish a "second order" spectral mapping theorem (Lemma [3|. 

In Section Isj we pursue a detailed analysis of accumulation points of the second order 
relative spectra on the real line. Our main contribution (Theorem [7| is a significant im- 
provement upon similar results previously found in (6] [7] • It allows calculation of rigourous 
convergence rates when the test subspaces are generated by a non-orthogonal basis. Con- 
crete applications include the important case of a finite element basis which was not 
covered by [3 Theorem 2.1]. Our present approach relies upon an homotopy argument 
which yields a precise control on the multiplicity of the second order spectral points. The 
argument is reminiscent of the method of proof of Goerisch Theorem; [151 122) . 

Theorem [Tj also determines the precise manner in which second order spectra encode 
information about the multiplicity of points in the spectrum of A. When an approximating 
space is "sufficiently close" to the eigenspace corresponding to an eigenvalue of finite 
multiplicity, a finite set of conjugate pairs in the second order spectrum becomes "isolated" 
and clusters near the eigenvalue. It turns out that the total multiplicity of these conjugate 
pairs exactly matches the multiplicity of the eigenvalue. This indicates that second order 
spectra detects in a reliable manner points in the discrete spectrum and their multiplicities, 
even under the variational collapse condition ([l]). In Section 14] we examine the practical 
validity of this statement on benchmark differential operators for subspaces generated by 
a basis of finite elements. 

The final section is aimed at finding the minimal region in the complex plane where the 
limit of second order spectra is allowed to accumulate. It turns out that, modulo a subset 
of topological dimension zero, this minimal region is completely determine by the essential 
spectrum of A. This gives an insight on the difficulties involving the problem of finding 
conditions on the test subspaces to guarantee convergence to the essential spectrum. 
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1.1. Notation. Below we will denote by Dom(A) the domain of A and by Spec(A) its 
spectrum. We decompose Spec(A) in the standard manner as the union of essential and 
discrete spectrum: 

Spec^^^(A) — {A G Spec(A) : 3x„ G Dom(A), ||a;„|| = 1, a;„ ^ 0, {A - A)a;„ -^ 0}, 
SpeCdi,(A) := Spec(^)\Spec,,,(A). 

The discrete spectrum is the set of isolated eigenvalues of finite multiplicity of A. 

Let (a, b) C R. Below D(a,6) will be the open disc in the complex plane with centre 
(a + 6)/2 and radius (6 — a)/2, and D[a,6] = D{a,b). We allow a — — oo or b = +cx) in 
the obvious way to denote half-planes or the whole of C. For b = a, D[a,a] — {a} and 
\D{a,a) = 0. 

Let A^ be a Hilbert space. Let fi C /C be an arbitrary subset and let Z3 C /C be a finite 
subset. We will denote 

disttc{t3,Q) := m&x ini \\u — v\\ic and distK:(w, fJ) = distK:({u}, fi). 

Here we include the possibility of AC = C and write dist(u, Q) — distc(u, f2). 
Let p G IN U {0}. We endow Dom(A'') with the graph inner product 

p 
{u, v)p ■- ^{A'u, A''v) Vu, V G Dom(A'') 

SO that (Dom(A^), {•, •)p) is a Hilbert space with the associated norm denoted by || • ||p. 
Below we will consider sequences of finite-dimensional subspaces {£„) = ('C„)„gN growing 
towards Dom(A'') with a given density property determined as follows; 

Ap = {(£„) C Dom(A^) : V/ G Dom(A''), distDom(AP)(/,/:n) ^ 0}. 

Note that if \\A\\ < oo, then Ao = Ap for any p G It^. If (jC-n) C Ap and £„ C jCn+i, then 
Uj^i^Cn is dense in the graph norm of A^. 

For a family of closed subsets n„ C C , the limit set of this family is defined as 

lim n„ = {z G C : 3zn G fin, Zn ^>- z} — {z £ C : dist(2:, fin) — > 0}. 
The limit set of a family of closed subsets is always closed. 

Remark 1. Below we will establish properties o/ lim„^oo Spec2(^,/In) for sequences 
(jC-n) G Ap. By similar arguments to those presented in this paper, one can establish 
analogies to all these properties for the "weak" limit set 



oo oo 



fl U Spec2(^,£„ 



n^l m^n 



if we impose that for all f G Dom(^'') exists a subsequence fn(k) £ C-nik) such that 

ll/n(fe)-/||p^O. 

Let bi, . . . ,bn be basis for a subspace C dH. Without further mention, we will identify 
the elements u £ C with a corresponding u G C" in the standard manner: 

n 

" = 5Z(w.^*)&i and u={{u,bl),...,{u,b*„)), 

fc=i 

where {fe*} is the basis conjugate to {6^}. If the bj are mutually orthogonal and ||6j|| = 1, 
then b* = bj. Below we will denote the orthogonal projection onto C hj P : H — > C. 
When referring to sequences of subspaces [Cn) we will use P„ instead. Note that {Cn) G Ao 
iff P„ — >■ / in the strong operator topology. 
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1.2. Toy models of spectral pollution in the Galerkin method. We now present 
a series of examples which will later illustrate some of the results below. See [21[ Theo- 
rem 2.1] and [25l Remark 2.5]. 

Example 1. Let H = span{e^}^i where e^ is an orthonormal set of vectors and let 
Cn = span{ef , . . . , e^_j^,e~}. Let /^ = ^(e;t ± e,T) and define 

in its maximal domain. Then 

Spec(P„A \ £„) = {0, ±1, . . . , ±(n - 1)} 
and so 

(2) lim Spec(P„A \ £„) \ Spec(y4) = {0} / 0. 

Note that the resolvent of A is compact. 

Example 2. If A is strongly indefinite and it has a compact resolvent, then there exists 
a sequence (Cn) G Ap for all p G IhJ such that 

(3) lim Spec(P„A \ £„) = R. 

Indeed, let Spec(y4) — {X^}mehi where the eigenvalues are repeated according to multi- 
plicity with A~^j^ < A" < and < A^ < A^^j. Let e^ be eigenvectors associated to 
A^ and assume that {e^} is an orthonormal set. Let {'ym.}jehi be an ordering of Q. For 
each n G It^ choose n < k € \f4 with A^ < 7j < A^ for 1 < j < n. Let 9j G (— tt, tv] be such 
that 

Ij = cos^(6'j)Afe+j + sin^(6lj)A^+j for 1 < j < n 
and define 

/„j = cos(6»j)efc+j. -f sin(6ij)e^+^^, 
then ([3| holds for the subspaces 

£„ = spanlej*", . . . ,ej,/„i,. . . , f„u} C Dom(A^). 

Example 3. Let "H and C„ be as in Example [l] Let 

/,t = sin(l/n)ej ± cos(l/n)eJ. 

For r > 0, define A = 5i;„g(.j fi''l/n )(/n 1 - J2„eM \fn){fn\ in its maximal domain. Then 

Spec^g^(A) = {-1} and Spec^^^{A) = {n''}„et~i. 

It is not difficult to see that now 

Spec(P„ A \ C„) = {-1, n sin(l/n)^ - cos(l/r^)^ 1, . . . , (ti - 1)"} 

where —1 is an eigenvalue of multiplicity n. Then, if r = 2, ([2| holds true. Note that the 
resolvent of A is not compact and A is now semi-bounded below. 

2. The second order spectra of a self-adjoint operator 

Definition 1 (See [21]). Given a subspace £ C Dom(j4), the second order spectrum of A 
relative to C is the set 

Spec2(A, C):={zeC : 3ue C\{0}, {{A - zl)u, {A - zl)v) = V?; G £}. 

If £ C Dom(yl^), then 

Spec2(^,£) = {zG C : 3u e C\{0}, {P{A - zfu,v) =0 Wv e £.} 

= {2 G C : G Spec(P(A - zf \ £.)}. 

Typically Spec2(A, £) contains non-real points. From the definition it is easy to see that 
z G Spec2(j4, C) if and only if I G Spec2(A, £). 
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2.1. Algebraic and geometric multiplicity. Let 

L — span{6j}"^i C Doiii(A) 

where the bj are linearly independent. Let B,L,M £ C"X" be matrices with entries given 
by 

(4) Bjk = {Abk, Abj), Ljk = {Abk, bj), M^k = (6fc, bj). 

With Q{z) = B — 2zL + z^M, it is readily seen that z £ Spec2(^, C) if and only if 

(5) 3m G C" \ {0} with Q{z)u = 0, 

that is, Spec2(A, yC) = Spec((5(2)). Evidently, Spec2(j4,ZI) is independent of the basis 
chosen, and since detM 7^ 0, it follows that Spec^iA, C) consists of at most 2n points. 
The quadratic eigenvalue problem ([5| may be solved via suitable linearisations, cf. [161 
Chapter 12] for details. For example, if 

then 2 G Spec2(yl, C) if and only if (T — zS)^ — for some v ^ 0. Equivalently, one can 
consider the matrix T where 

(6) ^'■^^''^^{-Ar'B 2i\/-% 
clearly Spec(T) — Spec{Q{z)) = Spec2{A,£). Consider also the non-singular matrices 

A straight forward calculation yields Q{z) (— /) — E{z){zl — T)F{z). The matrices 
Qiz) © {—!) and {zl — T) are said to be equivalent. Also, every matrix polynomial L{z) 
is equivalent to a diagonal matrix polynomial 

diag[ii(2), 12(2), . . . , ir{z), 0, . . . , O] 

where the diagonal entries ij{z) are polynomials of the form 11^(2; — Zjk) ''' with the 
property that ij{z) is divisible by ij-i{z). The factors (2 — Zjk)^'' are called elementary 
divisors and the Zjk are eigenvalues of L{z). The degrees of the elementary divisors 
associated to a particular eigenvalue zo are called the partial multiplicities of L{z) at 
2o; see [161 Theorem A. 6.1] for further details. For a linear matrix polynomial (such as 
T — zl) the degrees of the elementary divisors associated to a particular eigenvalue zq 
coincide with the sizes of the Jordan blocks associated to 20 as an eigenvalue of T; see 
[16| Theorem A. 6. 3]. The notion of multiplicity of an eigenvalue now follows from the 
fact that two n x n matrix polynomials are equivalent if and only if they have the same 
collection of elementary divisors; see [161 Theorem A. 6. 2]. 

Definition 2. The geometric and algebraic multiplicity of an element z £ Spec2(^,£) 
are defined to be the geometric and algebraic multiplicity of z as an eigenvalue ofT. 

Note that the definition of geometric and algebraic multiplicity is independent of the 
basis used to assemble the matrices B,L,M and T. Indeed, let B,L,M be assembled 
with respect a basis 61, . . . , 6„ and B, L, M be assembled with respect a basis Ci, . . . , c„, 
where C = span{6j}"^i = span{cj}"^i. Then, 

, 7 \ fN-' 0\/0 I \ [N 

■ ' -M~^B 2M-^l) \ N-^J y^Ar^B 2M'^LJ \0 N 

[cj,b*), and therefore T and T are equivalent. 
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2.2. The approximate spectral distance. Suppose that the given basis {bj}"^i oi £. 
is orthonormal so that A/ = J in Q. Let Q{z) = B — 2zL + z^I. For 2; £ C we define 

vec" \\v\\ 

Note that the right hand side is independent of the orthonormal basis chosen for C. If 
£ C Dom(4^), then 

, , . \\P{A-zfv\\ 
crA.c(z) = mm ,-— , 

vec \\v\\ 

and for z ^ Spec2(^,£) we have aA,ci^) ~ \\Qi^)^^\\~^ ■ Furthermore, z £ Spec2(j4,£) 
if and only if (Ja.c{z) = 0. The map aA,c '■ C — > [0, 00) is a Lipschitz subharmonic 
function; see jHi Lemma 1], [71 Lemma 4.1] and ^ Theorem 2]. Therefore Spec^iA, £) is 
completely characterised via the maximum principle. 

2.3. The spectrum and the second order spectra. Let 

Amin(^) = inf[Spec(yl)] and Amax(^) = sup[Spec(A)]. 
For all C C Dom(A), 

(7) SpeC2(yl, C) C D[An,in(A), An,ax(A)]; 

see [261 Theorem 3.1] and [28; Corollary 2.1]. Thus 

lim SpeC2(A,/:„) C D[An,in(A),An,ax(^)] 
n — ^00 

for any sequence {£n)- 

A growing interest in the second order relative spectrum and corresponding limit sets 
has been stimulated by the following property: if (a, b) n Spec(^) = 0, then 

(8) Spec2(A,£)nD(a,6) = V£ C Dom(yl); 
see [211 Theorem 5.2] or Lemma [I] below. Thus, 

(9) Spec(A) n [Re 2; - |Im 2:|,Rc 2 + |Im z|] / whenever 2 G Spec2(A,£); 

see [261 Corollary 4.2] and 21, Theorem 2.5]. Thus inclusions of points in the spectrum 
of A are achieved from Re 2 with a two-sided explicit residual given by |Im z\. 

In fact, the order of magnitude of the residue in the approximation of Spec (A) by 
projecting Spec2(^,£) into R can be improved to |Im 2p, if some information on the 
localisation of Spec(j4) is at hand. Indeed, if {a,b) n Spec(A) = {A} and z € Spec2(^,£) 
with 2 G D(a, b), then 



(10) 



|Im2p „ llmzl^ 

Re 2 - ' „' ,Re 2 + -i ^— 

— He 2 He z — a 



n Spec(yl) = {A}; 



see [TDl Corollary 2.6] and :28, Theorem 2.1 and Remark 2.2]. 

The following lemma is an improvement upon [261 Theorem 5.2]. It immediately implies 
property ([8|. 

Lemma 1. Let a,b £ R be such that (a, b) n Spec(yl) = 0. If z £ D{a,b), then 

(11) f^Axiz) > a{z) 
where 

(b — a) — \z — a\ — \z — b\'^ , , , ,,.9 

(12) Q(2)=Q(,,f,)(2)= ^ ' ^ '_ ' Ldist(2,{a,b})2 >0. 
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Proof. Without loss of generality assume that Im 2: > 0. The region 

F^ ■- {fi - z : n e Spec(A)} 

is contained in two sectors, one between the real line and the ray re"'*^ {r > 0), the other 
between the real line and re''''+^''\ where < 61 + 62 < n/2. We have cos(6li + 62) > 0, 
and by the cosine rule, 

./, n^ (b-af -\z-bf -\z-a\^ 
cos(6'i +62) = ^ ', ' „, L . 

Z\z — o\\z — a\ 

Now (Fz)^ — {w^ : w G F^} is contained in a sector between the rays re~*^ ^ and re'^ ^ . 
Therefore 

min{Re y : y £ e''-''''''^\F,f} = min{|a - z\\ \b - z\'^} cos(6»i + 6*2) = a{z). 

By virtue of the spectral theorem. 

Re e'-'^^-'^^'>{Q(z)v,v)cr^ = Re e'^'''-''''\{A - z)v, {A-z)v) 



'^"'-"'^ f{li-zfd{E^v,v) 



Re e 

> a(z)||u|| . 



D 



From ([8| it follows that 

D(a, b) n hm Spec2(>l, £n) = whenever (o, b) n Spec(A) = 

n — ^oc 

and 

Rn lim Spec2(A,£„) C Spec(yl). 

n — ^00 

We now consider two examples where 

(13) lim Spec2(A,£„) = Spec(y4). 

n—¥oo 

The first one has a compact resolvent while the second one has —1 in the essential spec- 
trum. 

Example 4. Let H, Cn and A be as Example [T] Then 

Spec.,{A,Cr^) = {±m,±l, . . . ,±(n - 1)}. 
Example 5. Let H, Cn and A be as Example [S] Then 

Spec2{A, £„) = {a„ ±i'y„,-l, 1, . . . , (n- 1)''} 

where a„ = n"" sin(l/n)^ — cos(l/n)^ and 7„ = (n"" + 1) sin(l/n) cos(l/n). The geometric 
multiplicity of the second order spectral point —1 is 71 — 1 and its algebraic multiplicity is 
2(n— 1). Asn— >-oo, the non-real point 

r -1 <r < 1 

ctn ± iju — >< — l±z r — 1 
(^ cxj r > 1. 

Hence (tTsl holds true for r 7^ 1. 



It is naturally expected that if {£n) approximate very fast a Weyl sequence for A £ 
Spec(j4), then 

Ae lim Spec2(j4, £„). 

n — ^00 

This intuition can be made rigorous through the following statement which appears to be 
of little practical use. 
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Proposition 2. Assume that \\A\\ < oo. Let A £ Spec(yl). Let (jCn) d H he such that 
n = AimCn- If there exists Xn G £n such that \\xn\\ ~ 1 and \\{A — A)a;n||'" ~ £n — >■ 0, 
then 

D(A - 5„, A + 5„) n Spec2(A £„) / 0, 5„ = 2(1 + \\A\\fen, 

for all n large enough. 

Proof. Let 7^ = T be as in ([6| for an orthonormal basis of £,„. If 

(14) \\{n-X)v\\c2.<e\\v\y,., 
then 

(15) D(A-(5,A + 5)nSpec2(A,£)/0 where 5 = {e\\Tn - \f"~'')^ ■ 

Indeed, suppose that the left side set i n (|15[ ) is empty. Then (7^ — A) G £2nx2n ^^ invertible 
and dist(A,Spec(7;j)) > S. Condition (|14[) iniphes IKT^i - A)~^|| > e. But HH Lemma 1] 

11^^^"^^ II- |dct(7;.-A)| < 55^ -^ ■ 

Hence ( |15| ) follows from the contradiction. 

Let w„ = a;„ ® Aa;„ so that \\vn\\ ~ Vl + A^. Then 

Il(r„-A)t^||c2„ ||P„(^-A)2x„|| ^ ||(A-A)|| ||(^-A)x„|| 



< 



Il!^llc2" V1 + A2 - V1 + A2 

Thus ([15]| holds for 

P-A||^ i_^ 

(1 + A2)4S 

The conclusion is a consequence of the identity \\Tn — A|| < (1 + || A||)^. D 

2.4. Mapping of second order spectra. 

Lemma 3. For a, b,c,d £ R, let f{w) — aw + b, g{w) — cw -\- d and F{w) = f{w)g{w)~^ . 
If ad / cb and g{Ay^ e Bin), then 

(16) zeSpec^iA,/:) ^=^ F{z) eSpec.2iFiA),9{A)C). 
Moreover, the algebraic and geometric multiplicities of z and F(z) are the same. 

Proof. Since g{A) preserves linear independence, a set of vectors {6j}"^i C £ is a basis for 
jC if and only if the corresponding set of vectors {g{A)bj}"^i is a basis for g{A)jC. Assume 
that {bj} is a basis for C and let bj — g{A)bj. It is readily seen that (bj)* = g{A)~^b*. 
Let u,v £ C and u = g{A)u and v = g{A)v. Then u = w G C" for any u £ C, obtaining 
the left side from {bj} and the right side from {bj}. Let _B, L, M be as in Q for A and 
{bj}. Let B, L, M be as in Q for F{A) and {bj}. Now, let z G Spec2{A,£.), and note 
that g{A)~^ G B{H) implies that g{z) / 0. We have 

{{B - 2zL + z'^M)u, v) = {{A- z)u, {A - z)v) 

= [ad - cb)-''g{z)\{F{A) - F{z))u, {F{A) - F{z))v) 

= [ad - cby'^g{zf{{B - 2F{z)L + F{zfM)u,v), 

since u, u G C" are arbitrary we deduce that 

(17) B-2zL + z^M = [ad - cby^g{zf{B - 2F{z)L + F(zfM). 

Thus F{z) G Spec2{F{A), g{A)C), and moreover, since the function {ad — cb)~'^ g{w)^ is 
analytic and non-zero at z, the algebraic and geometric multiplicities of z and F{z) are 
the same; see |16l Theorem A. 6. 6]. D 
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3. Accumulation of the second order spectrum and multiplicity 

3.1. Neighbourhoods of the discrete spectrum. Throughout this section we assume 
that 

(18) (a, b) n Spec(^) = {Ai < • • • < AJ C SpcCdi^(^). 
Let 

where £ > is small enough to ensure that rrij := Rank(lj) is equal to the multiplicity 
of Xj. Let m — X]^=i '^i be the total-multiplicity of the group of eigenvalues Ai, . . . , As. 
We will denote by £ the eigenspace associated to the group of eigenvalues in {a,b), that 
is, £ = Rangel(£,,6)(^). We fix an orthonormal basis B = \ujYj=i of £^, where the Uj are 
eigenvectors associated to the eigenvalues \j. 

Lemma 4. Let z € D(a, &). Let 

(19) 7 = 2\/5m(6 - af [(1 + max{|a|, |&|})^ + (& - a)^(2\/5s + 8)] . 
I] the subspace L C Dom(^^) is such that 

distDom(A2)(S,'C) < 5 < a(„,6)(2:)[dist(z,Spec^)] 7"\ 
then 

o-A,c{z) > ,^_ ,,^ {ai^a,b){z) [dist(2:, Spec A)] - ^7) > 0. 

Proof. Let ce IR\(a, 6). Define 

S 

A = A + ^{c- \j)lj r>om{A) =Dom(A). 
Then A is self-adjoint and Spec(yi) — {c} U Spec(^) \ {Ai, . . . , As}. Let 

s 

^W = -I](c- A,)(A, +c-2z)lj, 

then K{z) is a finite rank operator with range £ and 

{A-zf = {A-zLf ~K{z). 
Evidently, {A ~- z)"^ is invertible and for all u € ?^ we have 

\\v + (A-zy'^K{z)vf 

= l|lR\{a.6)t'f + U{a,h)V + {A - z)''^ {{A - zf -{A- z)^) 1(„,6) t; f 

= ||lR\(a.6)wf + \\{A- Z)-^{A - zfl(^^t,)vf 

s 

>Pi{zf\\v\\\ 

where /3i(2:) = min{l, dist(z, Spec(A))/|c — z\}. 

Let Q{z) — P(A — z)^ \c- For each eigenvector Uj we assign & Vj £ C such that 
||A^(wj — Wj)|| < (5 for p = 0, 1, 2. Note that in general Pmj /i'j,however H^Wj— UjH <5- 



|2 
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The hypothesis on 5 ensures that 

\\Q{z)vj - (i - zfu,\\ < \\P{A - zfiuj - v,)\\ + \c- zfWPuj - u, II 

< WiA- zf{Uj -Vj)]] +|c-2|^||PUj -Uj\\ 

s 

< \\[{A - zif + Y^{c - AO(A. + c - 2z)ti\{uj - vj)\\ + 5\ 

i = l 

< 5/32(2), 

where p2{z) = (1 + |2|)^ + s^,{z) + |c — z\^ and ^{z) — maxi<j<s |(c — \j){\j + c — 2z) 
Note that A satisfies the hypothesis of Lemmalll Let a{z) be given by (12 1. Then 

\\{A^zyu,-Q{z)-^Pu,\\ 

<\\{A- z)-\, - Q{zy\j\\ + \\Q{z)-^Pu, - Q{z)-\j\\ 

< II (c— zy Uj — Q(zy^Vj\\ + a{z)^^5 

<\\{c- z)-^v, - 0(2)^^,11 + {.a{z)-^ + |c - z\~-'^)5 

< a(2)-'||(c - zy^Q{z)v, - vj\\ + {a{z)-' + \c-z\-^)5 

< q(2)"'||(c - zy''Q{z)v, - uj\\ + (2q(2)"' + \c- z\-^)5 

< a{zy^\c- zy^\\Q{zy.j - (i- 2)^Uj|| + (2a(2)"^ + |c- 2p^)(5 

< Hz)S 

where /Ssiz) = a(2)-^|c - z\-^[l32{z) + 2|c - z\'^ + a{z)]. 
Thus 

||[(i-2)-^-Q(2)-lp]7^(2)H| 

s 

J=l 

s 

<J2\'^-xy2z-c-xMl{A-zy''-Qizy'p]iM\ 

TTl 

< ^^{z)\\v\\J2m- zy'' - Q(zy'p]u,\\ 
fc=i 

< 5m^(2)^3(2)||t;||. 
Hence, for any v £ C, 

\\p{A - zifv\\ = \\Q{zy + PK{zy\\ 
>\\Q{zy'r'\\v + Qizy'PK{zy\\ 

> a{z)[\\v +{A- 2)"^7f (2)1-11 - ||(i - zy^K{zy - Q(2)~'PA:(2)«||] 
>a(2)[/3i(2)'-<5m^(2)/33(2)]||«||. 

By fixing 

_ ( b+{b-a) = 2b-a if Re 2 > (a + 6)/2 
^ ^ \ a - (6 - a) = 2a - b if Re 2 < (a + b)/2, 

we get 



6 — a< |c — 2I < — ?^{b — a) 
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Thus 

dist(2:,Spcc(yl)) I ^ ^ r, f^tu ^2 

Pi[z) = , ^(z) < 2V5(o - a) , 

C — 2 

/32(z) < (l + max{|a|,|6|})^ + (b-a)^(2sV5 + 10/4) and a{z) < ^^ ~ J^' . 

Therefore 

(^A,c{z) > a(z)[j5i{zf - 5m^j.(z)l33{z)] 

a(z)[A\st{z, Spec(yl))]^ 5m2^/l>{b - af[l32{z) + 2\c - z]^ + a{z)] 
~ \c— zp \c — zP 

Q(2)[dist(2:, Spec(A))]^ - ^7 

|c-zP 

> ,|, ^ (Q(2)[dist(2,Spec(A))]^-^7). 

5|o — a|^ 

D 
Lemma 5. iei A G SpeCjj;j,(^) 6e an eigenvalue of multiplicity m and let 

dx = 2-'/^dist(A, Spec(A)\{A}). 
If£cC, then 

(20) Spec2(A,£)nD(A-dA,A + dA) = {A} 

and A as member of Spec2{A, C) has geometric multiplicity m and algebraic multtplicity 
2m. 

Proof. Assume that z G Spec2(y4,£) n D(A — c^a, A + dA) and z ^ \. Then there exists a 
normahsed u G •C\{0} such that 

(21) {{A~zI)u,{A~zI)w)=Q 



for all w G £. In particular {u,Uk) = for each ut G £. If we take w = n in |21|, we 
achieve 

\\{A-Re z)uf - jlmzl^lluf - 2i\lin z\{( A- Re z)u,u) = 0. 
Thus |Im z\ = II (^ - Re z)u|| / and {{A - Re z)u, u) = 0, so that 

\\{A - X)uf = |A - Re z\'^\\uf + ||(^ - Re z)uf + 2(A - Re z){{A - Re z)u, u) 



Since 



we get 



= |A - Re zf + |Irn zf < 2|Im z\\ 



dist A, Spec yl \ { A} < j— — ^ ] ^ {A- \)u , 

||(7-l(A-d;,,A + d;,))w|| 

, ., ^ s n ^ \/2|Im z\ 

M = ||(/-1(.-....)).|| < dist(A.Spec(A)\{A}) < '■ 



Statement (20 1 follows from the contradiction. 

Let {bi, . . . , bn} be a basis for £, and let B, L and M be the matrices Q associated to 
this basis. Let v be an arbitrary member of £. We have 

= {{A^ XI)uk, {A - XI)v))n = {{B - 2\L + \^M)u^,v)c^ 

so that (B - 2AL + \^M)uk = 0. Hence 



\t - Asl 



Uk 
\Uk 



/ A / / 

^B 2L I M 



X- 1-0. 
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Since det(S') 7^ 0, •^ (^ ), • • • , {r^ ) r s-rs eigenvectors of T corresponding to the eigen- 
value A. For k j^ j we have {Muk,Uj)c" = (ufe, Uj)ji = 0, so the vectors {ui, . . . , it^} are 
hnearly independent. There is a one to one correspondence between the eigenvectors of T 
associated to A and the eigenvectors oi A associated to A. Thus, the geometric multiphcity 
of A € Spec2{A,£) is equal to m. 

Now let us compute the algebraic multiplicity. Let v,w£ C" and u G Span{ufc}^i. If 



u 
Alt 



= [T-\I] 



I 



-A 



I 
I 



we have w = Au + u and —M ^ Bv + 2M ^ Lw — Xw = Am. Therefore 



-Bv + 2XLv - yMv = 2\Mu - 2Lu = 



so we deduce that 



V 



+ 



where v € Sp&n{uk}kLi- 



Suppose now that 



= [T~\I] 



We have w = Aw and -M ^Bv+2M ^Lw-Xw = u,theTeiore -Bv+2XLv-X^Mv = Mu. 
Thus 

{Mu, u) = -{Bv - 2XLv + X^Mv, u) = -{v, Bu - 2XLu + X^Mu) = 

so that u = 0. This ensures that spectral subspace associated with A as eigenvalue of T 
is given by 

SP'^^K Alt )'(^ )'•••'( A^)'(i 



Lemma 6. Let K. be a Hilbert space and A — {ui, . . . , Um} C AC be a set of orthogonal 
vectors with \\uj\\ic > 1 for 1 < j < ni. Let £ be an n-dimensional subspace of K. where 
n > m. Let Wj = Puj where P : AC — > £ is the orthogonal projection associated to £. 
Let Wj{t) = tWj + (1 — t)uj for t G [0, 1]. // ||wj — Mj|| < 1/^/rn for each 1 < j < m, then 
there exist vectors {wm+i, ■ ■ ■ , Wn} such that {wi(t), . . . , Wm{t), Wm+i, . . . , Wn} is linearly 
independent for all t G [0, 1]. 



Proof. If aiwi + • • • + OmWrn = where not all of the Oj 



0, we would have 



El 



= ||ai(Mi - Wi) H h a^niUm " Wm)||^ 

^\aj\\\uj - wj\\] < i^\aj\\ /m 

which contradicts Holder inequality. Hence necessarily {wi, . . . ,Wm} is linearly indepen- 
dent. Let {wm+i, ■ ■ ■ , Wn} be any completion of {wi, . . . , Wm} to a basis oi£. Let t G [0, 1]. 
Suppose now that aiWi(i) -I- • • • -I- amWm{t) + am+iWm+i + • • • + anWn — 0. Then 



= ^a,ii;,(i)+ Y. 



3=1 



J— m+1 



3=1 3=1 
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The two terms on the right-hand side are orthogonal and therefore each must vanish. As 
the set {wi, . . . ,Wn} is linearly independent, all aj = ensuring the conclusion of the 
lemma. D 

3.2. Main result. 

Theorem 7. Let a,b ^ Spec(yl) be such that condition | |18[ ) hold for a group of eigenvalues 
{Ai < ... < As} with corresponding multiplicities ruj. Let m = X^^^i "^j ^^ their total 
multiplicity. Let 

d = dist({a, 6}, Spec(A) \ {Aj}J=i)- 
Let K — d^/7 where 7 is defined by (19 1. Let 

< e < min < — —7- — — — , min -^ — —^ — - > where \q = a and \s+i = 6. 

y-m^i'^K^/^ o<j<s 2 J 

If C C Dom(A^) is such that 

distDom(A2)(fi,^) < t£^ 
for an orthonormal set of eigenfunctions B associated to {AjjJ^i, then 

s 

(22) Spec2(yl, C) n D(a, fe) C [j D{\j - e, Xj + e) 

and the total algebraic multiplicity of the points in Spec2(A, £,) n ID(Aj — e, Aj + e) is 2mj . 
Proof. Let a = a — d and b — b + d, then 

■ f / \ ■ / r \ b-a iff a + b\ 

mi a,- rAz) — mm < a,- zAz) : z — e H > 

zeD(a,b) ("''■'^ ^ ee(-n,7r]\ (''•'')^ ' 2 2 J 

To find the right hand side we may assume without loss of generality that a = —r and 
b — r, then for z — x + iy where x'^ + y'^ = r, we have 

(2r + 2df - \x + iy + r + d\'^ - \x + iy - r - d\'^ 
2\x + iy + r + d\\x + iy — r — d\ 

We assume that a; > as the case where i < can be treated analogously. A straightfor- 
ward calculation shows that 



, , iz; -^ zui — \x -t lu -f- I -r ti\ — x -t ly — r — u ,. , . , , ,,,2 

"(a.5) (^) = ^ oui.-.,/.^^lui.„_._^l ^di«t(a; + iy,{-r - d,r + d}f 



I . \ ,/ , n \ ■ \x + iy — r — d\ ,2 
mm a,~ t^(x + ly) — did + 2r) mm tt = d , 

0<:c<r ("-'''^ '^^ ^ ' 0<x<r | X + jj/ -|- r -f dj 

2,2 2,2 

y +x — r y +3: ^ r 

hence Lemma l4| applied to the interval (fi, 6) ensures (22 1. 

Let us prove the second part. In the notation of Lemma pi lot A — B, IC — Dom(A^) 
and define subspaces £t ~ span{iui(t), . . . , Wm{t),Wm+i, . • . , Wn} for t £ [0, 1]. Note that 
B is also an orthogonal set in Dom(A^) and ||Mj||Dom(A2) > 1- Since ne^ < m~^", we get 
dim£t = n for all t G [0, f]. Now, 

W"^] - Wj{t)\\oon,(A^) = i|l"i - ^jllDom(A2) = t\\{l - P)lij |1 Dom(A2) 

< distDon,(^2)(uj,£) for all t G [0, 1], 
thus dist]3o,-,-j/^2)(Z3,£t) < KE^, and by virtue of LemmaHlwe obtain 

s 

Spec2{A,£.t)nD{a,b) c[JD{Xj -e,\j +£) for all t G [0, 1]. 

Denote by Tit) the linearisation matrix defined as in ([6| for Ct. According to Lemma [s] 
the algebraic multiplicities of Aj G SpecT(l) is 27nj. If we let iTj{t) be the spectral 
projection associated to 

SpecT(t)n D(Aj -e,Aj -l-e). 
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then 

^.(ti)-^.(t2) = -^£_ _ (r(ii)-c)~'-(r(t2)-c)~'dc. 

Thus iTj{t) — )■ nj{l) as f — )■ 1. By virtue of [2UI Lemma 4.10, Section 1.4.6], we have 
Rank(7rj(0)) =Rank(7rj(l)). D 

An immediate consequence of Theorem IT] is that if a sequence of test subspaces £„ 
approximates a normahsed eigenfunction u associated to a simple eigenvalue A £ Spec(j4) 
at a given rate in the graph norm of A^ , 

then there exist a conjugate pair Zn,'z^ £ Spec2{A, £,„) such that |A — z„\ = 0((5^"). A 
simple example shows that this estimate is not necessarily optimal. 

Example 6. Let H — span{e,i}5JLo, ^ = '^''^\s-t){et\ and 

Cn = span{ei, . . . , e„_i, a„eo + /3„e„} 

where a^ + /3^ = 1 and /?„ — 5> 0. Then 

Spec2(yl,£„) = {l,...,n- 1,7„,7^} 

where 7„ ~ m/3„. On the other hand 

II A*'((a„ - l)eo + /?„e„)|| ~ n''/?^ P = 0, 1, 2. 

Thus |7„| ~ n/3„ while distoojjj(^2)(S, £n) ~ ri^Pn- Moreover, note that in general 
distDojn(^2)(S, C„) might diverge as n -> cxd and still we might be able to recover |7„| — >■ 0; 
take for instance /3„ — 1/n^". 

An application of LemmalSlyields convergence to eigenvalues under the weaker assump- 
tion jCn e Al . 

Corollary 8. Let —oo < a < b < oo be such that (a, b) n SpeCj,j,g(A) — 0. If (Cn) G Ai, 
then 

(23) lim Spec2(A,£„) n D(a,6) = SpeCji^(yl) n {a,b). 

n — ^oo 

Proof. Choose e > such that a + e < min{A G Spec(yl) n (a, b)} and let fi £ {a,a + e). 
With g(w) — w ~ fi, the operator g{A)~^ is bounded and g{w) satisfies the hypothesis of 
Lemmals] Let Mn = g{A)C„- Then (Mn) G Aq. Indeed, any u £ % can be expressed 
as 11 = g{A)v for v G Dom(A), and since {Cn) G Ai, we can find «„ G Cn such that 
g{A)vn — > (7(j4)w. By virtue of Theorem ItI applied to g{A)~^ , 

lim Spec2(5(^)-\M„)nD(p(fe)-\3(a + e)-i) = SpeCdi,(<?(4)-i)n(p(b)-\3(a + e)-i). 

The fact that g{ui) is a conformal mapping, Lemma|3]and the spectral mapping theorem 
ensure the desired conclusion. D 

Example 7. Let A be the operator of examples^ and [5] By virtue of ||8| and CoroUaryls] 

lim Spec2(A,/:„) C Spec(yl) U (-1 + iR) 

n— J-oo 

for any {C„) G Ai. 

Corollary 9. If A has compact resolvent and (Cn) G Ai, then 

lim Spec2{A, Cn) = Spec(yl) 

n—^oo 

When A is not semi-bounded, this statement is in stark contrast to Example [2] 
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4. Numerical applications 

Calculating second order spectra leads to enclosures for discrete points of Spec(^). 
This can be achieved by combining ||9| with Theorem u\ The latter yields a priori upper 
bounds for the length of the enclosure in terms of bounds for the distance from the test 
space, £, to an orthonormal basis of eigenfunctions B. In practice, these upper bounds 
are found from interpolation estimates for bases of C. 

4.1. On multiplicity and approximation. Let A G Spec(A) be an isolated eigenvalue of 
multiplicity m, so that ( |22[ ) ensures an upper bound on the estimation of A by m conjugate 
pairs of Spec2(^, jC). li£. approximates well a basis of only I < m eigenfunctions and rather 
poorly the remaining m — l elements of B, then it would be expected that only I conjugate 
pairs of Spec2{A, C) will be close to A while the remaining m — l will lie at a substantial 
distance from this eigenvalue. We illustrates this locking effect on the multiplicity in a 
simple numerical experiment. 

Let A — —A — —d^ — dy subject to Dirichlet boundary conditions on f7 = [0, vr]^ C R^. 
Then 

Spec(A) = {f +k'^ ■.j,kehl} 

and a family of eigenfunctions is given by Ujk{x, y) — sm{jx) sm{ky). Some eigenvalues of 
A are simple and some are multiple. In particular the ground eigenvalue 2 = 1 + 1 is simple 
and 1 + fc^ for fc = 2, 3, 4, 5 are double. The contrast between the two eigenfunctions Uki 
and uik will increase as k increases. The former will be highly oscillatory in the a:-direction 
while the later will be so in the y-direction. If £ captures well oscillations in only one 
direction, we should expect one conjugate pair of Spec2(A, £) to be close to 1 + fc'^ and 
the other to be not so close to this eigenvalue. 

In order to implement a finite element scheme for the computation of the second order 
spectra of the Dirichlet Laplacian, the condition C C Dom{A) prescribes the corresponding 
basis to be at least C^-conforming. We let C — C{s) be generated by a basis of Argyris 
elements on given triangulations of fl. Contrast between the residues in the interpolation 
of uik and Uki is achieved by considering triangulation that are stretched either in the x 
or in the y direction. See Figure [l] right. 
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Figure 1. Left: change in the length of the enclosure predicted by (|9| 
for the eigenvalues A = 1 + fc'^ of the 2-dimensional Dirichlet Laplacian 
on [0,7r]^. The test subspaces C, are generated by Argyris elements on 
a uniform mesh of 240 triangles with a prescribed numbers of segments 
intersecting the boundary. Right: Corresponding mesh for the first 6 
iterates of the left graph. 
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In this experiment we generate triangulations with a fixed total number of 240 elements, 
resulting from diagonally bisecting a decomposition of fi as the union of 120 rectangles 
of equal ratio. We consider mesh with s = 3, 4, 5, 6, 8, 10, 12, 15, 20, 24, 30 and 40 
elements of equal size on the lower edge [0, n] x {0} C dQ. We then compute the pairs 
zik{s), zik{s), Zki{s), Zkiis) e Spec2{A, C{s)) which are closer to the eigenvalues 1 + fc^ 
than to any other point in Spec(A). By virtue of (|9|, we know that Re zik{s) — |Im zik{s)\ 
and Re 2a;i(s) — jIm zik{s)\ are lower bounds for this eigenvalue and Re 2ifc(s) + jlm 2ifc(s)j 
and Re Zki(s) + |Im zik(s)\ are corresponding upper bounds. 

Remark 2. Even though we have the same number of elements in each of the above mesh, 
in general they do not have the same amount of elements intersecting dQ. Then typically 
dim£(s) 7^ dim£(i) for s ^ t, although these numbers do not differ substantially. The 
precise dimension of the test spaces are as follows: 



s 


3 


4 


5 


6 


8 


10 


dim£(s) 


2774 


2648 


2578 


2536 


2494 


2480 



On the left side of Figurelllwe have depicted the residues 2|Im 2ifc(s)| and 2|Im Zki{s)\ 
for each one of the eigenvalues 1 + fc'^ in the vertical axis, versus s in the horizontal 
axis on a semi-log scale. The graph suggests that the order of approximation for all the 
eigenvalues changes at least two orders of magnitude as s varies. The minimal residue in 
the approximation of the ground eigenvalue is achieved when s = 10 and s = 12. This 
corresponds to low contrast in the basis of £(s). When the eigenvalue is multiple, however, 
the minimal residue is achieved by increasing the contrast in the basis. As this contrast 
increases, one conjugate pair will get closer to the real axis while the other will move away 
from it. The greater the k is, the greater contrast is needed to achieve a minimal residue 
and the further away the conjugate pairs travel from each other. 

This experiment suggests a natural extension for Theorem [T] If only I < m mem- 
bers of B are close to C, then only I conjugate pairs on Spec^iA, C) will be close to the 
corresponding eigenvalue. 

4.2. Optimality of convergence to eigenvalues. We saw in Example[6]that the upper 
bound established in Theorem [T] is sub-optimal. We now examine this assertion from a 
practical perspective. 

Let A = —d^ + V acting onl-L — L'^{Q,-k) where V is a smooth real- valued bounded 
potential. Let 

Dom(A) ^{ue H'\Q, tt) : u(0) = u{-k) = 0}, 
so that A defines a self-adjoint operator semi-bounded. Note that A G Spec(A) = 
SpeCjjig;. {A) if and only if A solves the Sturm-Liouville eigenvalue problem An — Xu subject 
to homogeneous Dirichlet boundary conditions at and tt. 

Let H be an equidistant partition of [0,7r] into n subintervals /; — [xi-i,xi] of length 
h — tt/u = xi — xi-i. Let 

£{h,k,r) = Vh{k,r,S) ^ {v £ 0^0,71) : v \i,e Pr{Ii), 1 < I < n, v{0) = v{n) ^ 0} 

be the finite element space generated by C'^-conforming elements of order r subject to 
Dirichlet boundary conditions at and tt; [T^ . An implementation of standard interpo- 
lation error estimates for finite elements combined with Theorem [T] ensures the following. 

Lemma 10. Let A, T-i and £{h, k, r) be as in the previous paragraphs. Let 

Spec(A) = {Ai < Aa < ...}. 

Let Oj = |Aj -I- |Aj_i and bj = jAj -I- |Aj+i, where Ao = — oo. For all r > k > 3, there 
exist a constant c > 0, dependant on j, k and r, but independent of h, such that 

Spec2{A,£.{h,k,r)) n \D{aj,bj) = {zhkr,zKi^} C D(Aj - ch^,Xj +ch'~^) 
for all h > sufficiently small. 
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Proof. Use the well-known estimate 

where Vh € C{h, k, r) is the finite element interpolant of u G C'^ (1 //^^^(O, tt); [121 Theo- 
rem 3.1.6]. Note that all eigenvalues of A are simple and its eigenfunctions are C°° . D 

Therefore each individual eigenvalue Xj is approximated by second order spectral points 

T--3 „ 

at a rate 0{h ^ ) for test subspaces generated by a basis of C -conforming finite elements 
of order r > 4. Due to the high regularity required on the approximating basis, this 
results is only of limited practical use. In fact, only A; > 1 (and r > 3) is required for 
C{h, k, r) C Dom(yl). Simple numerical experiments confirm that the exponent predicted 
by Lemma [To] is not optimal, see Figure [2] 




Table: p such that |Im Zhir{j)\ ~ h'' 



X 


r = 3 


r = 4 


r = 5 


1 


1.9979 


2.9900 


3.9849 


4 


1.9983 


2.9878 


3.9191 


9 


2.0023 


2.9838 


3.9235 


16 


2.0099 


2.9764 


3.9152 


25 


1.9962 


2.9673 


3.8871 



Element size h 

Figure 2. Figure corresponding to Example [8] Left: Logiog plot of 
2|Im Zhirij)\ versus element sizes. Right: Slopes of linear interpolations 
of these graphs. 



Example 8. Suppose that 1/ = so that A = — 9j. In this experiment we fix an 
equidistant partition S and let C = £{h, l,r) be the space of Hermite elements of order 
r = 3, 4, 5 satisfying Dirichlet boundary conditions in and n. We then find the conjugate 
pairs {zhirij), Zhir{j)} € Spec2(^,-C) which are close to Xj = j^ € Spec(j4) for j = 
1,2,3,4,5. On the left of Figure [2] we have depicted 2|Im Zhir{j)\ versus small values of 
h. On the right side we have tabulated the slopes of linear fittings of the lines found. 

This experiment suggests that the order of approximation of A € Spec^ji^j. (A) from its 
closest z £ Spec2{A, jC) is 

(24) |z-Al~distDom{A)(S,/:). 

Indeed, interpolation by Hermite elements of order r has an iif '^-error proportional to ft"""^ . 
The same convergence rate is confirmed by Example [6] 

Example 9. Let V{x) = 2cos(2a::) be the Mathieu potential. Let C be as in Example [s] 
The exponents p reported in Figure [s] further confirm (24 1. 



Remark 3. The error in the estimation of the eigenvalues of a one- dimensional elliptic 
problem of order 2p by the Galerkin method using Hermite elements of order r is pro- 
portional to h'^^^^^~^' ; [571 Theorem 6.1]. If A is a second order differential operator, 
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the quadratic eigenvalue problem ([5| gives rise to a non-self-adjoint fourth order prob- 
lem which is to be solved by a projection-type method. Thus (241 is consistent with this 



estimate if we take into account the improved enclosure (10 1, see Section 4-3 

4.3. Improved accuracy. Estimate (fTol can be combined with (|9| to provide improved 
a posteriori enclosures for A G Specjis;.(^). The key idea is to estimate an upper bound 
a> V = sup{A £ Spec(j4) : A < A}, a lower bound b <v = inf{A G Spec(A) : A > A} and 
2 G Spec2(j4, C) such that 

(Re z - |lm 2:|,Rc z + |lm z\) n Spec(A) = {A}. 

Here ^, and v can be elements of the discrete or the essential spectrum of A. The one- 
sided bounds a and b can be found from Spec2(A,/I) or by analytical means. If 6 — a is 
sufficiently large and |Im z\ is sufficiently small, (10 1 improves upon ([9}. We illustrate 
this approach a practical settings. 

Example 10. Let A = -dl + 2cos(2a;), H = L^(0,7r) and C = £{h,l,r) for r = 3,4,5 
be as in Example [9] In Figure [S] we have computed inclusions for the first five eigenvalues 
of A by directly employing ([9| and by the technique described in the previous paragraph. 
In the latter case, we have found Oj — a from the computed upper bound for A-,-i < Xj 
{ai — — oo) using (|9|. Similarly for bj — b. These calculations can be compared with 
those in [2 Tables 1,2,3]. See also JS) §7.4]. 
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Figure 3. Enclosures for \j in the case oi A = —d^ + 2cos(22:) by 
direct application of ([9| and by the improved ( 10 1. Here C = C{-,l,r). 



We now consider an example from solid state physics to illustrate the case where p and 
v are not in the discrete spectrum. 

Example 11. Let A = —d^ + cos(a;) — e^^ acting on L^(— co,co) where Dom(A) = 
H^{—oo, oo). Then A is a semi-bounded operator but now SpeCggg(^) consists of an infinite 
number of non-intersecting bands, separated by gaps, determined by the periodic part of 
the potential. The endpoints of these bands can be found analytically. They correspond to 
the so called Mathieu characteristic values. The addition of a fast decaying perturbation 
gives rise to a non-trivial discrete spectrum. Non-degenerate isolated eigenvalues can 
appear below the bottom of the essential spectrum or in the gaps between bands. 

In Figure |4] we report on computation of inclusions for the first three eigenvalues in 
SpeCjj;j,j,(A): Ai < minSpeCj.g3(j4), A2 in the first gap of the essential spectrum and A3 in 
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the second gap. No other eigenvalue is to be found in any of these gaps. Here C, is the 
space of Hermite elements of order 3 on a mesh of 20i segments of equal size ft = 0.1 in 
[— /, Z] subject to Dirichlet boundary conditions at ±L For the improved enclosure, a and 
h are approximations of the endpoints of the gaps where the Aj lie. 
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Figure 4. Enclosures for the first three eigenvalues of the perturbed 
periodic Schrodinger operator of Example |11| The numerical values of 
a and h are found in [l] (see also |101 Table 1]). 

Since the eigenvectors of A decay exponentially fast a& x -^ ±cx), the members of 
SpeCjj;^^. (A) are close to eigenvalues of the regular Sturm-Liouville problem Au = \u 
subject to u{—l) = u{l) = for sufficiently large I < cxa. The numerical method considered 
in this example does not distinguish between the I = 00 and the large / < 00 eigenvalue 
problem. For instance, the inclusion found for Ai in Figure [4] is also an inclusion for the 
Dirichlet ground eigenvalue of A \ L^{—25, 25). 

In the case of A2 and A3, which lie in gaps of SpeC;,j,j, (A) , they should be close to high 
energy eigenvalues of the finite interval problem. Indeed, for the parameters considered in 
FigureUI the inclusion for A2 is also an inclusion for the 17"" eigenvalue of A f L^(— 50, 50). 
Similarly that for A3 is an inclusion for the 65"^ eigenvalue of A fL^(— 100, 100). 

Remark 4. The error in the Galerkin approximation of the j-th eigenvalue of a regular 
Sturm-Liouville problem via Hermite elements of order 3 on an uniform mesh of element 
size h is known to he 0{j^h^); .27> §6.2-(34)]. Evidently the latter is only an upper bound, 
there might be high energy eigenvalues which are approximated accurately: such as the 
17*'' eigenvalue of A\L'^ {~50, 50) or the 65"" eigenvalue 0/ A |"L^(-100, 100) m the above 
example. In practice it is not easy to take advantage of this observation as the Galerkin 
method also produces a considerable amount of spurious eigenvalues in low-energy regions 
of the spectrum. These correspond to approximations of singular Weyl seguences of points 
in the essential spectrum of A. 

5. Accumulation points outside the real line 

Since R is second countable and the essential spectrum of A is closed, there always 
exists a family of open intervals {aj,bj) C R such that R \ SpeCj,j,j,(yl) — U^i(aj, 6j). 
Throughout this section we will be repeatedly referring to the following two yl-dependent 
regions of the complex plane: 

00 
B = Ba = [J ID(aj,6j) and A = Aa = D[inf Spec^,,^, supSpec^^^^] \ B. 

For any given set SI C C and e > we will denote the open e -neighbourhood of SI by 

[Q]e = {zG C : dist(2,!r2) <e}. 

Lemma 11. LetN = A/" C B be compact and such thatAfr]Spec{A) = 0. Let {£n) G A2. 
There exists sj^ > and N > such that 

o"A,£„(2) > sj\f \lz £ N and n> N. 
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Proof. B = if and only if Spec(A) — R and in this case there is nothing to proof, so 
we assume B 7^ 0. Since A/" is compact and every covering of a compact set has a finite 
sub-covering there exists a finite family J" such that A/" C Uj^j^D{aj,bj). Then we can 
decompose Af = VJj^tNj where each Nj = A/j are compact and such that A/j C ID(aj, hj). 
Since A/'nSpec(j4) = 0, there exists e > such that Aj C D(aj-|-e, bj—e). By construction 
each interval {aj, bj) can only intersect Spec(yl) at a finite number of isolated eigenvalues 
of finite multiplicity. Therefore the conclusion follows from lemmas [T] and |4] applied in 
{oj + e, bj — e) using analogous arguments as in the proof of TheoremTTJ D 

Combining this lemma with a similar argument as in the proof of Corollary [S] we 
immediately achieve the following theorem which generalises Corollary |9] 

Theorem 12. For any {£n) G Ai, 

Spec<ji3,(A) C lim Spec2(yl,£„) C A U SpeCdi,e(^). 

n—i-oc 

Let us now examine statements complementary to this result. We begin with two 
auxiliary lemmas. 

Lemma 13. Let c £ Spec j,^^ (A). For any given m £ \t\ and 5 > Q there exists {xi\YL\ C 
l(c_5,c+5)(^) such that 

(i) ||a;i|| = 1 for all I — 1, ... ,m. 

(ii) {APxi,A''xi) = for all I / [ and p,q = 0,1 

(iii) Axi — cxi + xi where \\xi\\ < 5 for any I — 1, . . . ,m. 

Proof. If c is an isolated point of the spectrum, then c is an eigenvalue of infinite multiplic- 
ity and the proof is elementary. Otherwise, there exists ci ^ c such that c; £ Spec(^)\{c}. 
By substituting ci by a sub-sequence if necessary, we can assume that |c; — c| = S{1) < S/2 
for 0<S{1 + 1) < (5(0 such that 



/ S{1) ^S{1)\ 



By picking xi £lf _a(o ,Hi)^\{A) for I — 1, . . . ,m with ||x;|| = 1 the desired conclusion 
is guaranteed. D 

For the triplet ci < C2 < C3 denote 

A(ci,C2,C3) = D[C1,C3]\ (D(C1,C2)UD(C2,C3)). 

Lemma 14. Let |c|<l,0<S<7r and \c\ < a < 1. Let z — ae'* and 



If z e A(-l,c, 1), then l3± > 0. 

Proof. The proof involves straightforward trigonometric arguments. D 

Lemma 15. iei 2 £ A(ci, C2, C3). For S > let xi,X2,X3 (^Doin{A) be such that 

(i) \M\ = l 

(ii) {APxk,A''x^^) =Qfork^k and p,q^ 0,1 
(iii) Axk = CkXk + Xk where \\xk\\ < S. 
There exist at G R independent of S such that if y ~ X]fe=i f^kXk, the polynomial 

{Ay,Ay)-2\{Ay,y)+\'{y,y) 
has roots A-|- = 2 + 0{5) and A_ = 2 + 0(5). 
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Proof. The cases where Ck — cj. for k ^ k are easy, so we assume that ci < C2 < C3. 
Moreover, without loss of generahty we can assume that ci = — 1, C2 = c where jcj < 1 
and C3 = 1. Then z and z axe roots of the polynomial 

a — 2Xa cos 6 + X 

for z = ae'*. If y = Y^l^^ o-kXk we get 

(y,j/) = Q? + 0:2 + ^3 

(-42/, y) = -a? + CQ2 + 03 + 0((5) 
{Ay,Ay)^al + ^al+al + 0{5). 
Finally note that the solution of the system 

"i + "2 + Q3 = 1 
— Qi + ca2 + 03 = a cos 9 

2,22, 2 2 

Qfi + C a2 + Ct-i, ^ a 



is easily found to be ai = \//?+ , 0^2 = \ ^_°2 a-nd 03 = ^/^^ where /3± are the expressions 
found in Lemma [141 D 

Theorem 16. Let N = N ^ k he compact. For all e > there exists C^ C Dom(j4) such 
that 

(25) JVc[Spec2{A,C,)]^c[JVhe. 

Proof. Let A/" and e > be fixed. Since A/" is compact, there exists {zj}"^i C Af such that 
Af C [(2j)"=i]_^ ,^ C [A/']e/2- The proof wiU be completed if £.. C Dom(A) is such that 

(26) [{^.}?=i],/4 ^ [Spec2{A,C,)l C [{^.}"=i],/, ■ 

Below we will choose the parameter S > small enough. 

Let {cjk}kZl'^jZl C Spec^gg(A) be such that Zj G A(cij, C2j,C3j) for all j = 1, . . . , n. By 
virtue of Lemma 
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there exists a family of vectors {xjk}k^i'''j^i ^ Dom(A) such that 
(i) ||a;fcj|| ~ 1 for all fc = 1, 2, 3 and j = 1, . . . ,n 
(ii) {A''xkj,A''x-^-.) = for all {k,j) / (fc, j) and p, g = 0, 1 
(iii) Axkj = CkjXkj + Xkj where \\xkj \\ < 5 for any fc = 1, 2, 3 and j = 1, . . . , n. 
Let akj = Q-k be as in Lemma [T5| for z = Zj. Choosing (5 > small enough and defining 
3 
Vi = ^ akjXkj and C^ = span{yi, . . . , y„}. 



yields ( 26 1 and hence ( 25 1 . D 



This result has two straightforward consequences. Given any compact subset N = 
A/" C A, there exists {£n} C Dom(yl) such that M = lim„^oo Spec2(A, £„). Evidently 
(£„) might fall outside Aq in general. On the other hand, we can always find {J~.n) £ Ai 
such that 

lim Spec(A, £„) = A U SpeCdig^(A). 
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